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The crystalline or liquid character of the downward cusp states in A'^-electron parabolic quan- 
tum dots (QD's) at high magnetic fields is investigated using conditional probability distributions 
obtained from exact diagonalization. These states are of crystalline character for fractional fillings 
covering both low and high values, unlike the liquid Jastrow-Laughlin wave functions, but in re- 
markable agreement with the rotating- Wigner-molecule ones [Phys. Rev. B 66, 115315 (2002)]. 
The crystalline arrangement consists of concentric polygonal rings that rotate independently of each 
other, with the electrons on each ring rotating coherently. We show that the rotation stabilizes the 
Wigner molecule relative to the static one defined by the broken-symmetry unrestricted-Hartree- 
Fock solution. We discuss the non-rigid behavior of the rotating Wigner molecule and pertinent 
features of the excitation spectrum, including the occurrence of a gap between the ground and first- 
excited states that underlies the incompressibility of the system. This leads us to conjecture that 
the rotating crystal (and not the static one) remains the relevant ground state for low fractional 
fiUings even at the thermodynamic limit. 
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I. INTRODUCTION 

The excitation energy spectrum of a two-dimensional 
A^-electron semiconductor quantum dot (QD), plotted as 
a function of angular momentum for a given high mag- 
netic field [B) , exhibits downward cusp a^i^i?i^i^ at certain 
magic angular momenta (L,„), corresponding to states 
with enhanced stability. For a given value of B, one of 
these Lm 's corresponds to the global minimum of the en- 
ergy, that is to the ground state (the ground-state value 
of Lm is denoted as Lg^)- Varying the magnetic field 
causes the ground state and its angular momentum Lgg 
to change. We note that due to their enhanced stability, 
only cusp states can become ground states. Underlying 
these properties is the inherent incompressibility of the 
cusp states in response to an external magnetic field. As a 
result, the cusp states have been long recognizedii2i^i^*Sii 
as the finite- A'^ precursors of the fractional quantum Hall 
states in extended systems. In particular, the fractional 
fillings I' (defined in the thermodynamic limit) are related 
to the magic angular momenta of the finite- A^ system as 
follows!^ 



(Henceforth, we will drop the subscript m, unless neces- 
sary). 

In the literature of the fractional quantum Hall effect 
(FQHE), ever since the celebrated paper— by Laughlin 
in 1983, the cusp states have been considered to be the 
antithesis of the Wigner crystal and to be described accu- 
rately by liquid-like wave functions, such as the Jastrow- 
Laughlin^'iS (JL) and composite-fermion»ii^ (CP) ones. 
This view, however, has been recently challenged^^ by 
the explicit derivation of trial wave functions for the 
cusp states that are associated with a rotating Wigner 



(or electron) molecule, RWM. As we discussed^ earlier, 
the parameter-free RWM wave functionSfi^ which are by 
construction crystalline in character, promise to provide 
a simpler, but yet improved and more consistent descrip- 
tion of the properties of the cusp states, in particular for 
high angular momenta (corresponding to low fractional 
fiUings). 

Issues pertaining to the liquid or crystalline charac- 
ter of the cusp states are significant in both the fields of 
QD's and the FQHE. Since the many-body wave func- 
tions in the lowest Landau level (high B) obtained from 
exact diagonalization (EXD), the RWM wave functions, 
and the CF/JL ones have good angular momenta™ L > 
Lq = N{N — l)/2, their electron densities are circularly 
symmetric. Therefore investigation of the crystalline or 
liquid character of these states requires examination of 
the conditional probability distributions (CPD's, i.e., the 
fully anisotropic pair correlation functions). These cal- 
culations were performed here under high magnetic field 
conditions for QD's (in a disk geometry^^) with A^ = 6—9 
electrons, and for an extensive range of angular mo- 
menta. This allowed us to conclude that in all instances 
examined here (corrresponding to 0.467 > u > 0.111) 
the cusp states exhibit an unmistakably crystalline char- 
acter, in contrast to the long held perception in the 
FQHE literature, with the RWM yielding superior agree- 
ment with the exact-diagonalization results^ Further- 
more, the RWM states are found to be energetically 
stabilized (i.e., exhibit gain in correlation energy) with 
respect to the corresponding static (symmetry-broken) 
Wigner molecules, from which the multideterminantal 
RWM wavcfunctions are obtained through an angular- 
momentum projection. We will present arguments that 
allow us to conjecture that the stabilization energy of the 
cusp states in high B remains nonvanishing even in the 
thermodynamic limit. 
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In the beginning of sect. II, we display the Hamilto- 
nian of the system under consideration and define the 
conditional probability distributions. Subsequently, in 
the same section, we present our main results pertaining 
to the structural properties of the CPD's. Possible im- 
provementes of the RWM wave functions are discussed in 
sect. 11(c). In sect. Ill, we recapitulate the essential as- 
pects of the two-step method of symmetry breaking and 
symmetry restoration, calculate stabilization energies for 
the RWM, and discuss pertinent features of its excitation 
spectrum (in particular, the occurrence of a gap between 
the ground state and the first excited state that is not 
a mere consequence of finite size; the appearance of this 
gap underlies the incompressibility of the system). A 
discussion pertaining to implications for the thermody- 
namic limit is presented in sect. IV. Finally, a summary 
is given in sect. V. 

II. CONDITIONAL PROBABILITY 
DISTRIBUTIONS 

We are interested in wave functions which are exact 
solutions (or good approximations to them) of the two- 
dimensional many-body problem defined by the Hamil- 
tonian 

JV n ^ * N N y 

2—1 i—1 i—1 jyi 

(2) 

which describes N electrons (interacting via a Coulomb 
repulsion) confined in a parabolic potential of frequency 
loq and subjected to a perpendicular magnetic field _B, 
whose vector potential is given in the symmetric gauge 
by 

A(r) = iBxr=i(-By,Ba;,0); (3) 

m* is the effective electron mass, k is the dielectric con- 
stant of the semiconductor material, and = jr^ — rj|. 
For sufficiently high magnetic field values (i.e., in the 
FQHE regime), the electrons are fully spin-polarized and 
the Zeeman term (not shown here) does not need to be 
considered. 

In the B oo limit, the external confinement luq can 
be neglected, and H can be restricted to operate in the 
lowest Landau level (LLL), reducing to the form: 

,r'i^c \ - \ ^ e 

i—l j>z 

where u)c = eB/{m*c) is the cyclotron frequency. 

For finite N, the solutions to the Schrodinger equa- 
tions corresponding to the Hamiltonians Q and Q must 
have a good angular mome ntum L. As described by 
us in detail in Refs. 17 (see also sect. Ill be- 

low) , these solutions can be well approximated by a two- 
step method of symmetry breaking at the unrestricted 



Hartree-Fock (UHF) level and of subsequent symmetry 
restoration via post-Hartree-Fock projection techniques. 
As elaborated in our earlier work^^StiSiii the two-step 
method describes the formation and properties of rotat- 
ing Wigner molecules in QD's. 

As indicated in the Introduction, probing of structural 
characteristics in many-body wave functions with good 
angular momentum L requires the use of the conditional 
probability distributions defined by 

N N 

P(r, r„) ^{^l\J2Y. ^(^^ - '•)'^(^^- - ro)|<f L>/(<f L|<i>L). 

(5) 

Here ^^(ri, r2, r^r) denotes the many-body wave 
function under consideration. In this paper, we calculate 
the CPD's for three types of many-body wave functions 
defined in the lowest Landau level: (i) the analytic rotat- 
ing Wigner molecule wave function, $RWM ^-g^g g^jg^ ggj,_ 
tion III); (ii) the wave function ^f^^ obtained through 
exact diagonalization in the lowest Landau level; and (iii) 
the Jastrow/Laughlin functions ^j^- 

P(r,ro) is proportional to the conditional probability 
of finding an electron at r under the condition that a 
second electron is located at Tq. This quantity positions 
the observer on the moving (intrinsic) frame of reference 
specified by the collective (coherent) rotations that are 
associated with the good angular momenta of the cusp 
states. 

A. Crystallinity in lower fractions l/9<2^<l/5 

The CPD's for cusp states corresponding to a lower 
filling factor than ly — 1/5, calculated for = 6 electrons 
with L = 135 {v = 1/9, left column) and for A^ = 6 with 
L = 105 (j/ = 1/7, right column) are displayed in Fig. 1. 
Fig. 2 displays the CPD's for the cusp states with A = 6 
electrons and L ~ 75 {v = 1/5, left column) and N = 7 
and L — 105 {i^ — 1/5, right column). In both figures, the 
top row depicts the RWM case. The EXD case is given 
by the middle row, while the CF cases (which reduce to 
the JL wave functions for these fractions) are given by 
the bottom row. 

There are three principal conclusions that can be 
drawn from an inspection of Figs. 1 and 2: 

(I) The character of the EXD states is unmistakably 
crystalline with the EXD CPD's exhibiting a well devel- 
oped molecular polygonal configuration [(1,5) for A = 6 
and (1,6) for A = 7, with one electron at the center], in 
agreement with the explicitly crystalline RWM case. 

(II) For all the examined instances (covering the frac- 
tional fillings 1/9, 1/7, and 1/5), the JL wave func- 
tions fail to capture the intrinsic crystallinity of the EXD 
states. In contrast, they represent "liquid" states in 
agreement with an analysis that goes back to the origi- 
nal papersSiiS by Laughlin. In particular, Ref . Hol investi- 
gated the character of the JL states through the use of a 
pair correlation function [usually denoted by g{R)] that 
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FIG. 1: Conditional probability distributions at high B for 
A*' = 6 electrons with L — 135 {i/ = 1/9, left column) and 
L = 105 {f = 1/7, right column). Top row: RWM case. Mid- 
dle row: The case of exact diagonalization. Bottom row: The 
Jastrow-Laughlin case. It is apparent that the exact diagonal- 
ization and RWM wave functions have a pronouned crystalline 
character, corresponding to the (1,5) polygonal configuration 
of the rotating Wigner molecule. In contrast, the Jastrow- 
Laughlin wave functions fail to capture this crystalline char- 
acter, exhibiting a rather "liquid" character. The observa- 
tion point (identified by a solid dot) was placed at the maxi- 
mum of the outer ring of the radial electron density^ of the 
EXD wave function, namely at tq = 7.318/_b for L — 135 and 
ro = 6AA21b for L = 105. Here, Ib = {hc/eBY^'^. The EXD 
Coulomb interaction energies (in the lowest Landau level) are 
1.6305 and 1.8533 b^/kIb for L — 135 and L — 105, respec- 
tively. The errors relative to the corresponding EXD energies 
and the overlaps of the trial functions with the EXD ones are: 
(I) For L = 135, RWM: 0.34%, 0.860; JL: 0.50%, 0.665. (II) 
For L = 105, RWM: 0.48%, 0.850; JL: 0.46%, 0.710. 






FIG. 2: Conditional probability distributions at high B for 
A'' = 6 electrons and L = 75 {v = 1/5, left column) and 
for N = 7 electrons and L = 105 (again v — 1/5, right 
column). Top row: RWM case. Middle row: The case of 
exact diagonalization. Bottom row: The Jastrow-Laughlin 
case. The exact diagonalization and RWM wave functions 
have a pronouned crystalline character, corresponding to the 
(1,5) polygonal configuration of the RWM for N = 6, and 
to the (1,6) polygonal configuration for N = 7. In contrast, 
the Jastrow-Laughlin wave functions exhibit a characteristic 
liquid profile that depends smoothly on the number A'^ of elec- 
trons. The observation point (identified by a solid dot) is lo- 
cated at ro = 5.431/s for N = 6 and L = 75 and ro — 5.883/s 
for N = 7 and L = 105. The EXD Coulomb interaction en- 
ergies (lowest Landau level) are 2.2018 and 2.9144 c^/kIb 
for iV = 6, L = 75 and iV = 7, i = 105, respectively. The 
errors relative to the corresponding EXD energies and the 
overlaps of the trial functions with the EXD ones are: (I) For 
iV = 6,L = 75, RWM: 0.85%, 0.817; JL: 0.32%, 0.837. (II) 
For iV = 7, L = 105, RWM: 0.59%, 0.842; JL: 0.55%, 0.754. 



determines the probability of finding another electron at 
the absolute relative distance R = |r — ro| from the ob- 
servation point Tq. Our anisotropic CPD of Eq. Q is 
of course more general (and more difficult to calculate) 
than the g{R) function of Ref. (lO|. However, both our 
P(r, ro) (for = 6 and N = 7 electrons) and the g{R) 
(for N = 1000 electrons, and for = 1/3 and ly = 1/5) 
in Ref. ^3 reveal a similar characteristic liquid-like and 
short-range-order behavior for the JL states, eloquently 
described in Ref. (see p. 249 and p. 251). Indeed, we 
remark that only the first-neighbor electrons on the outer 
rings can be distinguished as separate localized electrons 
in our CPD plots of the JL functions [see Fig. 1 and Fig. 
2]. 

(Ill) For a finite number of electrons, pronounced crys- 
tallinity of the EXD states occurs already at the rather 
high 1/ = 1/5 value (see Fig. 2). This finding is par- 
ticularly interesting in light of expectationsi^ (based on 
comparisongSiiSiiS between the JL states and the static 
bulk Wigner crystal) that a liquid-to-crystal phase tran- 
sition may take place only at lower fillings with < 1/7. 



B. Crystallinity in higher fractions 1/5 < < 1 

Following the conclusion that the crystalline charac- 
ter of the cusp states in QD's is already well devel- 
oped for fractional fillings with the unexpected high value 
of — 1/5, a natural question arises concerning the 
presence or absence of crystallinity in cusp states cor- 
responding to higher fractional fillings, i.e., states with 
1/5 < V < 1/3, and even with 1/3 < v < 1. To an- 
swer this question, we show in Fig. 3 the CPD's for the 
RWM (left column) and EXD (right column) wave func- 
tions for the case of A^ = 7 electrons and for L = 69 
(1/5 <v = 7/23 < 1/3), i = 57 (1/3 <u = 7/19 < 1), 
and L = 45 (1/3 < v = 7/15 < 1). Unlike the long 
held perceptions in the FQHE literature (which were re- 
asserted in two recent publicationsiS) the CPD's in Fig. 
3 demonstrate that the character of the cusp states with 
high fractional fillings is not necessarily "liquid-like" . In- 
deed, these high-i^ cusp states exhibit a well developed 
crystallinity associated with the (1,6) polygonal configu- 
ration of the RWM, appropriate for N = 7 electrons. 

Of interest also is the case oi v — 1/3. Indeed, for this 
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FIG. 3: Conditional probability distributions at high B for 
N = 7 electrons and L ^ 69 {ly = 7/23 = 0.304 > 1/5, 
top row), L ^ 57 {I > u = 7/19 = 0.368 > 1/3, middle 
row), and L = 45 (1 > 1/ = 7/15 = 0.467 > 1/3, bottom 
row). RWM case: Left column. The case of exact diago- 
nalization is depicted in the right column. Even for these 
low magic angular momenta (high fractional fillings), both 
the exact-diagonalization and RWM wave functions have a 
pronouned crystalline character [corresponding to the (1,6) 
polygonal configuration of the RWM for N = 7 electrons]. 
The observation point (identified by a solid dot) is located 
at ro = 4.7521b for L = 69, ro = 4.2781b for L = 57, and 
ro = 3.7761b for L = 45. 



N=7, L=63 (v=1/3) 




FIG. 4: CPD's at high B for iV = 7 and L = 63 (i^ = 1/3). 
Top: RWM case; Middle: EXD case; Bottom: JL case. Un- 
like the JL CPD (which is liquid), the CPD's for the exact- 
diagonalization and RWM wave functions exhibit a well devel- 
oped crystalline character [corresponding to the (1,6) polyg- 
onal configuration of the RWM for N = 7 electrons]. The 
observation point (identified by a solid dot) is located at 
ro = 4.5681b. 



fractional filling, the liquid JL function is expected to 
provide the best approximation, due to very high overlaps 
(better than 0.99) with the exact wave function. ■^^■^^ In 



Fig. 4, we display the CPD's for N = 7 and L = 63 
{i^ = 1/3), and for the three cases of RWM, EXD, and JL 
wave functions. Again, even in this most favorable case, 
the CPD of the JL function disagrees with the EXD one, 
which exhibits clearly a (1,6) crystalline configuration in 
agreement with the RWM CPD. 

Similar crystalline correlations at higher fractions were 
also found for QD's of a different size, e.g., with = 6, 

= 8, and A^ = 9 electrons. As illustrative examples for 
these additional sizes, we display in Fig. 5 the CPD's for 
AT = 8 audi = 91 (1/5 4/13 < 1/3) and for A^ = 9 

and L = 101 ( 1/3 < 1/ = 36/101 < 1). Again, the CPD's 
(both for the RWM and the EXD wave functions) exhibit 
a well developed crystalline character in accordance with 
the (1,7) and (2,7) polygonal configurations of the RWM, 
appropriate for A^ = 8 and A^ = 9 electrons, respectively. 

The case of A^ = 9 is of particuler significance. In- 
deed it represents the smallest number of electrons with 
a non-trivial concentric-ring arrangement, i.e., the inner 
ring has more than one electrons. As the two CPD's 
(reflecting the choice of taking the observation point [ro 
in Eq. (O] on the outer or the inner ring) for A^ = 9 
reveal, the polygonal electron rings rotate independently 
of each other. Thus, e.g., to an observer located on the 
inner ring, the outer ring will appear as uniform, and 
vice versa. The fact that both the RWM and exact wave 
functions share this property of independently rotating 
rings is a testament to the ability of the RWM theory to 
capture the essential physics of QD's in high B. 



C. Improvements of the RWM wave functions 

It is of importance to note here that the favorable com- 
parison between the crystalline structure of the RWM 
and that of the exact wave functions, in contrast to the 
liquid-like character of the JL functions, persists even 
for cases where the latter is found to have the advan- 
tage (over the RWM) concerning total energies and wave 
function overlaps. As examples, we refer to the case of 
v = 1/5 discussed in the caption of Fig. 2 (see in par- 
ticular the errors in the energies for the RWM and JL 
functions relative to the exact energies given at the end 
of the caption). 

Close inspection of the humps in the CPD's obtained 
from the RWM and through exact diagonalization reveals 
that the RWM tends to somewhat overestimate the de- 
gree of crystallinity, i.e., the RWM humps are narrower 
and higher (this tendency diminishes for larger values of 
L). Nevertheless, the degree of overall agreement be- 
tween the exact results and those obtained through the 
parameter-free RWM wave functions is rather remark- 
able. Moreover, the high level of agreement between the 
RWM and exact results extends to other properties. This 
includes the zeroes (often called vortices) of the many- 
body wave functions. Indeed, as recently shown in Ref. 
123 . the exact wave functions (in contrast to the JL ones) 
have simple zeroes whose topology is in agreement with 
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FIG. 5: Additional CPD's at high B. RWM resuhs: Left col- 
umn. Results from exact diagonalization are depicted on the 
right column. Top row: N = 8 electrons and L — 91 (1/5 < 
ly = 4/13 = 0.308 < 1/3). Two bottom rows: = 9 elec- 
trons and L = 101 (1/3 <u = 36/101 = 0.356 < 1, see text 
for explanation). Even for these low magic angular momenta 
(high fractional fillings), both the exact-diagonalization and 
RWM wave functions have a pronouned crystalline character 
[corresponding to the (1,7) and (2,7) polygonal configurations 
of the RWM for A'' = 8 and 9 electrons]. The observation 
point (identified by a solid dot) is located at ro = 5.105Zfl for 
A' = 8,L = 91, and ro = 5.218/s (outer) and ro = 1.662Zfl 
(inner) for iV = 9,L = 101. 



that of the simple zeroes of the RWM functions. 

The above suggest that the RWM wave functions can 
form the nucleus for constructing a whole class of rotat- 
ing crystalline functions with added variational freedom, 
which will yield further quantitative energetic and struc- 
tural improvements. For example the RWM functions 
could be used as the basis for constructing variational 
wave functions in diffusion^'^ and variational^^ quantum 
Monte-Carlo studies. For a most recent investigation 
along these lines, see Ref. where our RWM function 
is augmented by a Jastrow prefactor with an exponent 
that is treated variationally. We remark, however, that 
the variational wave function employed in Ref. |25| has 
multiple zeroes due to the Jastrow factor, in disagree- 
ment with the exact diagonalization results. 



this point requires generation of global ground states out 
of the cusp states, achieved through inclusionSiii of an 
external parabolic confinement (of frequency ujo). In the 
two-step method, the static WM is first described by an 
unrestricted Hartree-Fock (UHF) determinant that vio- 
lates the circular symmetry^ Subsequently, the rotation 
of the WM is described by a post-Hartree-Fock step of 
restoration of the broken circular symmetry via projec- 
tion techniquesii^ We note that, in the limit N — > oo, 
the static WM of the UHF develops to the extended two- 
dimensional Wigner crystal^i and its more sophisticated 
variants. 

In general, the localized broken symmetry orbitals of 
the HF determinant are determined numerically via a 
selfconsistent solution of the UHF equations. Since we 
focus here on the case of high B, we can approximate 
the UHF orbitals (first step of our procedure) by (pa- 
rameter free) displaced Gaussian functions; namely, for 
an electron localized at (Zj), we use the orbital 
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cxp 



2A2 



(6) 



with A — yjT%/m*Q.; = ^/oj^+lo^JA, where uJc = 
eB/{m*c) is the cyclotron frequency. We have used 
complex numbers to represent the position variables, so 
that z = X + iy, Zj = Xj + iYj. The phase guaran- 
tees gauge invariance in the presence of a perpendicular 
magnetic field and is given in the symmetric gauge by 
(p{z,Zj;B) = {xYj - yXj)/2ll, with Ib = ^ficjeB. We 
only consider the case of fully polarized electrons, which 
is appropriate at high B. 

We take the Zj 's to coincide with the equilibrium po- 
sitions [forming a structure of r concentric regular poly- 
gons denoted as (rti, rt2, «r)] of = X]q=i clas- 
sical point charges inside an external parabolic confine- 
ment of frequency Wq- Then we proceed to construct the 
UHF determinant ^^^^[z] out of the orbitals u{zi, .Zi)'s, 
i — 1,...,N. Correlated many-body states with good 
total angular momenta L can be extractec&i& from the 
UHF determinant using projection operators, i.e.. 
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c*™^(7i,...,7.)exp[^*^7,^9j , (7) 



III. RESTORATION OF CIRCULAR 
SYMMETRY 

A. Correlated many-body wave functions 

Our two-step method for deriving the RWM wave func- 
tions is anchored in the distinction'^^ between a static and 
a rotating Wigner molecule, with the rotation stabilizing 
the latter relative to the former. Further elaboration on 



where L = J2q=i-^q 4"^^^ [7] is the original UHF 
determinant with all the orbitals of the qth ring rotated 
(collectively, i.e., coherently) by the same azimuthal an- 
gle 7q, and each ring is rotated independently of each 
other. Note that Eq. lO can be written as a product of 
projection operators acting on the original UHF deter- 
minant [i.e., ^''^^^(71 — 0,...,7r = 0), see Eqs. (6) and 
(7) in Ref. Setting A = IbV2 restricts the orbital in 
Eq. ijHJl to lie entirely in the lowest Landau level, and al- 
lows for the derivation of the analytic RWM functionsi^ 
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FIG. 6: Stabilization energies ASff" for iV = 6 (dashed 
curve) and N = 1 (solid curve) fully polarized electrons in 
a parabolic QD as a function of B. The troughs associated 
with the major fractional fillings (1/3, 1/5, and 1/7) and the 
corresponding ground-state angular momenta [see Eq. Q] are 
indicated with arrows. We have extended the calculations up 
to B = 120 T (not shown), and verified that ASIJ'" remains 
negative while its absolute value vanishes as _B ^ oo. The 
choice of parameters is; Tiuq = 3 meV (parabolic confine- 
ment), m* = 0.067me (electron effective mass), and k = 12.9 
(dielectric constant). 



We stress that while the initial trial wave function of the 
UHF equations consists of a single determinant, the pro- 
jected wave function is a linear superposition of many 
determinants, as can be explicitly seen from the analytic 
forms of the RWM functions in Refs. (S). 



B. Stabilization energy 

In the case of finite B (requiring the inclusion of con- 
finement, i.e., Wo 7^ 0), the projected energy correspond- 
ing to a symmetry-restored RWM state with angular mo- 
mentum L is given [in the case of a single (0, TV) or a 
(l,iV- 1) ring] ^ ^ 
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7i(7)e'^^d7, (8) 



with h{j) = (*UHF(o)|iy|,I,UHF(^)^ ^j^j ^ 

(^fUHF(0)|*UHF(^)^^ where *uhF(^) jg original UHF 
determinant with all the orbitals rotated (collectively) by 
the same azimuthal angle 7. H is the many-body Hamil- 
tonian in Eq. The UHF energies are simply given 

by £;uHF = /i(0)/n(0). 

We note that, unlike the UHF ground state (describ- 
ing a static Wigner molecule) which does not have good 
angular momentum, the ground states of the RWM ex- 
hibit good angular momenta (labeled afore- 
mentioned) that coincide with magic ones [we denote the 
ground-state energy of the RWM as -Ef^j = -BpRj(igs)]- 
Note that in Fig. 6 the ground-state magic angular mo- 
menta obey 



with fc = 0, 1, 2, ... Such sequences, having as a period 
the number of eletrons on the crystalline polygonal ring 
[5 and 6 for the (1,5) and (1,6) RWM's corresponding to 
iV = 6 and = 7] , reflect directly the collective rotation 
and incompressibilty of the RWM (see Sect. III.C below). 



The stabilization energy, Ai?g'™ = E\ 
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E 



UHF 



Xg, = N{N - l)/2 + k{N - 1), 



(9) 



the rotating WM relative to the static one [namely the 
fact that iJf'^j < i?uHF, see Fig. 6] is a purely quantum 
effect. This energy gain, Ai?||"", demonstrated here for 
= 6 and 7 electrons, is in fact a general property of 
states projected out of trial functions with broken sym- 
metry. This is due to an "energy gain" theoremS^ stat- 
ing that at least one of the projected states (i.e., the 
ground state) has an energy lower than that of the orig- 
inal broken-symmetry trial function (e.g., the UHF de- 
terminant considered above). 



C. Excitation gap 

The oscillations of the stabilization energy [Fig. (6)] 
reflect the oscillatory behavior of the energy of the pro- 
jected RWM states, as well as of the exact ones, since the 
mean-field energy E'uhf varies smoothly with B (see Ref. 
[T^ . Underlying the oscillatory behavior of the ground- 
state energies is a fundamental property of the spectrum 
of the system, namely, the appearance of special gaps due 
to the enhanced stability of the cusp states. Indeed, for a 
given magnetic field, both the ground state (specified by 
Lgs), as well as the first excited state (specified by Li), 
are magic [withi^ Li = Lgs ± (A^ - 1) for TV = 6 - 8]. 
As an example of this behavior, we display in Fig. 7 the 
low part of the EXD excitation spectrum for a QD with 
N ^7, B = 18.8 T, hujQ = 3 meV, and k = 12.9. The 
states L — 99, L — 105 (the ground state), and L = 111 
are demonstrated indeed to be cusp states of enhanced 
stability (all three states are well separated from the rest 
of the excited states). 

For the given magnetic field B = 18.8 T, the first ex- 
cited state corresponds to L = 99. However, as B in- 
creases, we found that the state with L = 111 diminishes 
in energy relative to that with L ~ 99, becoming itself 
the first excited state, and eventually (as B is increased 
further) replacing L — 105 as the ground state. This 
sequence of changes occurs for all ground states (with 
magic angular momenta) accessed through variation of 
the magnetic field. This results in the behavior of the 
excitation gap, A = — E^^, shown in Fig. 8 [the cal- 
culations here were performed with the RWM projected 
energies of Eq. (jS))]. The gap in Fig. 8 separates states 
with similar internal structure, i.e., they exhibit the same 
polygonal configuration as revealed through the CPD 
analysis for N — 7 electrons (see Figs. 2 — 4). The inter- 
nal structures of higher excited states differ from that of 
the ground state, with the disparity increasing with the 
excitation energy. 

The incompressibility of the cusp states (which, as dis- 
cussed above, correspond to magic angular momenta) is 



FIG. 7: Low-energy part of the spectrum of the parabolic 
QD whose parameters are the same as those in Fig. 6, cal- 
culated as a function of the angular momentum L through 
exact diagonalization for N = 7 electrons with a magnetic 
field B — 18.8 T. We show here the spectrum in the interval 
95 < L < 115 (in the neighborhood of = 1/5). The magic 
angular momentum values corresponding to cusp states are 
marked (99, 105, and 111), and they are seen to be separated 
from the rest of the spectrum. For the given value of B, the 
global energy minimum (ground state) occurs for Lgs — 105, 
and the gap A to the first excited state (L = 99) is indicated. 
The lowest energies for the different L's in the plotted range 
are connected by a dashed line, as a guide to the eye. The 
zero of energy corresponds to 7hQ, where D, = {luq -1-0)^/4)^^^ 
and LOc ~ eB/{m*c). The horizontal arrow denotes the en- 
ergy of the Laughlin quasihole at L = 112. It is seen that the 
Laughlin quasihole is not the lowest excited state. 



connected directly to the appearance of the gap (as dis- 
cussed in the context of the FQHE in Ref. 9). The dis- 
cussion presented here about the nature of the excitation 
spectrum (and in particular the existence of a gap A) 
allows us to comment on the influence of impurities and 
disorder on the properties of the quantum dot. Natu- 
rally, we focus on the regime of small or moderate disor- 
der, since a high degree of disorder (or strong impurities) 
will destroy both the gaps in the spectrum, as well as 
the coherent (collective) nature of the rotating Wigner 
molecule. The effect of disorder (or impurities) depends 
on the size of the gap. For sufficiently weak disorder, both 
the excitation gap A and the coherence of the magic an- 
gular momentum states maintain, namely, the states sep- 
arated by the gap experience only local disorder-induced 
perturbations (i.e., they broaden) and they remain con- 
ducting (see Sect. IV below). Obviously, for cases of a 
vanishing gap (i.e., between the fractional fillings, see ar- 
rows in Fig. 8), even weak disorder can induce a global 
change in the character of the (perturbed) wave functions 
by strongly mixing the degenerate Lgg and Li cusp states 





10 11 12 13 14 

B (T) 

FIG. 8: Top: RWM projected energies [see Eq. 0] calculated 
as a function of the magnetic field B for N — 7 electrons in a 
parabolic QD with the same parameters as those used in Fig. 
6 and Fig. 7. Each of the parabola-like curves (made partly 
of a solid and partly of a dashed line) corresponds to the 
marked value of the angular momentum - i.e., for the range 
of magnetic fields shown here, Lgs = 57 (7/19), 63 (1/3), 
69 (7/23), 75 (7/25), and 81 (7/27), with the corresponding 
value of the fractional filling v = N{N — l)/(2Lgs) given in 
parentheses. The solid lines denote the ground-state energies, 
and the dashed lines give the values of the first-excited-state 
energies. Note that the gap between the ground and the first 
excited state, A = — _E^°, oscillates as a function of B. 
The arrows denote the values of B for which the gap vanishes, 
occuring between the fractional fillings. The zero of energy 
corresponds to 7hQ + _E^f, where Q ^ {luq + 0)^/4)^^^ [with 
ujc ~ eB/{m*c)] and E^l is the classical energy of the static 
Wigner molecule (see Ref. IT^ . Bottom: The gap A plotted 
versus B. 

(and often additional nearby cusp states depending on 
the magnitude of B), and this can lead to a state with a 
broken-symmetry electron density having characteristics 
of a pinned Wigner crystallitei^ 



IV. DISCUSSION: IMPLICATIONS FOR THE 
THERMODYNAMIC LIMIT 

While our focus here is on the behavior of trial and 
exact wave functions in (finite) QD's in high magnetic 
fields, it is natural to inquire about possible implications 
of our findings to FQHE systems in the thermodynamic 
limit. 

We recall that appropriate trial wave functions for 
clean FQHE systems possess a good angular momentum 
L> Lq, a. property shared by both the CF/JL and RWM 
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functions i^i^iifliii We also recall the previous finding^- 
that for large fractional fillings v > vq, the liquid-like 
(and circularly uniform) CF/JL function is in the ther- 
modynamic limit energetically favored comparedSiiSiiSiiS 
to the broken-symmetry static Wigner crystal (which has 
no good angular momentum); for v < vq, the static 
Wigner crystal becomes lower in energy. This finding 
was enabled by the form of the JL functions, which fa- 
cilitated computations of total energies as a function of 
size for sufficiently large N (e.g., N = 1000). 

A main finding of this paper is that the exact- 
numerical-diagonalization wave functions of small sys- 
tems (N < lOj are crystalline in character for both low 
and high fractional fillings. This finding contradicts ear- 
lier suggestion82i2iiSii& that, for high i^'s, small systems 
are accurately described by the liquid-like JL wave func- 
tions and their descendants, e.g., the composite-fermion 
ones. Of course, for the same high i/'s, our small-size 
results cannot exclude the possibility that the CPD's of 
the exact solution may exhibit with increasing N a tran- 
sition from crystalline to liquid character, in agreement 
with the JL function. However, at the moment the ex- 
istence of such a transition remains an open theoretical 
subject. 

For the low fractions, the RWM theory raises still an- 
other line of inquiry. Due to the specific form of the 
RWM wave functions, computational limitations'^^ pre- 
vent us at present from making extrapolations of total 
energies at a given v as N ^ oo. Nevertheless, from the 
general theory of projection operators, one can conclude 
that the RWM energies exhibit a different trend com- 
pared to the JL ones, whose energies were found^iiSiiSiiS 
to be higher than the static Wigner crystal. Indeed the 
rotating- Wigner-molecule wave functions remain lower in 
energy than the corresponding static crystalline state for 
all values of N and z/, even in the thermodynamic limit. 
This is due to the fact that the aforementioned energy- 
gain theorem.^^ (see sect. Ill) applies for any number of 
electrons N and for all values of the magnetic field B. 
Naturally, the RWM wave functions will be physically 
relevant compared to those of the broken-symmetry crys- 
tal at the thermodynamic limit if the energy gain does 
not vanish when oo; otherwise, one needs to con- 

sider the posssibility that the static crystal is the relevant 
physical picture. 

The discussion in the above paragraph may be reca- 
pitulated by the following question: which state is the 
relevant one in the thermodynamic limit [N — > cx)) - 
the broken-symmetry one (i.e., the static Wigner crys- 
tal) or the symmetry restored (i.e., rotating crystal) 
state? This question, in the context of bulk broken- 
symmetry systems, has been addressed in the early work 
of AndersonpS, who concluded that the broken-symmetry 
state (here the UHF static crystalline solution) can be 
safely taken as the effective ground state. In arriving at 
this conclusion Anderson invoked the concept of (gener- 
alized) rigidity. As a concrete example, one would ex- 
pect a crystal to behave like a macroscopic body, whose 



Hamiltonian is that of a heavy rigid rotor with a low- 
energy excitation spectrum I? j^J , the moment of in- 
ertia J being of order A^ (macroscopically large when 
A'^ — > cxd). The low-energy excitation spectrum of this 
heavy rigid rotor above the ground-state (L — 0) is es- 
sentially gapless (i.e., continuous). Thus although the 
formal ground state posseses continuous rotational sym- 
metry (i.e., L = 0), "there is a manifold of other states, 
degenerate in the N ^ oo limit, which can be recombined 
to give a very stable wave packet with essentially the 
nature"— of the broken-symmetry state (i.e., the static 
Wigner crystal in our case). As a consequence of the 
"macroscopic heaviness" as A^ ^ oo, one has: (I) The 
energy gain due to symmetry restoration (i.e., the stabi- 
lization energy A£'||''") vanishes as A^ ^ oo, and (II) The 
relaxation of the system from the wave packet state (i.e., 
the static Wigner crystal) to the symmetrized one (i.e., 
the rotating crystal) becomes exceedingly long. This pic- 
ture underlies the aforementioned conclusion that in the 
thermodynamic limit the broken-symmetry state may be 
used as the effective ground state. 

Consequently, in the following we will focus on is- 
sues pertaining to the "rigidity" of the rotating Wigner 
molecule in high magnetic fields. In particular, using our 
projection method and exact diagonalization, we demon- 
strated explicitly in previous publications^^''^'^ that the 
rigid-rotor picture applies to an A^-electron QD only 
when B — 0. In contrast, in the presence of a high mag- 
netic field, we found^^ that the electrons in the QD do not 
exhibit global rigidity and therefore cannot be modeled 
as a macroscopic rotating crystal. Instead, a more ap- 
propriate model is that of a highly non-rigid rotor whose 
moment of inertia depends strongly on the value of the 
angular momentum L. 

The non-rigid rotor at high B has several unique prop- 
erties: (I) The ground state has angular momentum 
Lgs > 0; (II) While the rotating electron molecule does 
not exhibit global rigidity, it possesses azimuthal rigidity 
(i.e., all electrons on a given ring rotate coherently), with 
the rings, however, rotating independently of each other. 
Furthermore, the radii of the rings vary for different val- 
ues of L, unlike the case of a rigid rotor (see for example 
the locations of ro in Fig. 1 and Fig. 3 for different L 
values); (III) The excitation spectra do not vary as L^; 
instead they consist of terms that vary as aL + b/\n^ 
[for the precise values of the constants a and b in the 
case of (0, A^) or (1, A^ — 1) electron molecules, see Ref. 
[T^ ; (IV) The angular momentum values are given by the 
magic values [see Ref. L = + X]g=i ^g"-?' where 
(ni, 712, n^) is the polygonal ring arrangement of the 
static Wigner molecule (with rig the number of electrons 
on the qth ring) and ki < k2 < ... < kr are nonnegative 
integers. These magic L's are associated with the cusp 
states which exhibit a relative energy gain with respect to 
neighboring excitations. Thus the low-energy excitation 
spectrum of the non-rigid rotor is not dense and exhibits 
gaps due to the occurrence of the magic (cusp) states 
(see section III.C). Furthermore, these gaps are reflected 
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in the oscillatory behavior of A£'||"" (see, e.g., Fig. 6) as 
a function of B (or ly) . 

As N increases, more polygonal rings are successively 
added, and since the polygonal rings rotate indepen- 
dently of each other (see, e.g., the case of iV = 9 in 
Fig. 5), we expect that the non-rigid- rotor picture re- 
mains valid even as iV — > oo. As a result, it is plausi- 
ble to conjecture the following properties at high B in 
the thermodynamic limit: (I) the oscillatory character of 
Ai?!!"" will maintain, yielding nonvanishing stabilization 
energies at the fractional fillings [see Eq. (QJ], and (II) 
the low-energy excitation spectra of the system will still 
exhibit gaps in the neighborhood of the magic angular 
momenta (see Fig. 8). Of course, these conjectures need 
to be further supported through numerical calculations 
for large N. Nevertheless, the above discussion indicates 
that the question of which state is physically relevant for 
low fractions in the thermodynamic limit at high B - i.e., 
the broken-symmetry static crystal or the symmetrized 
rotating crystal ~ remains open, and cannot be answered 
solely following the path of Ref. IIE 

The rotating Wigner crystal has properties characteris- 
tic of FQHE states, i.e., it is incompressible (connected to 
the presence of an excitation gap) and carries a current^Si 
(while the broken-symmetry static crystal is insulating). 
Thus, we may conjecture that a transition at lower frac- 
tional fillings from a conducting state with good circular 
symmetry to an insulating Wigner crystal cannot occur 
spontaneously for clean systems. Therefore, it should be 
possible to observe FQHE-type behavior at low fractional 
fillings in a clean system - a prediction that could ex- 
plain the observations of Ref. |33, where FQHE behavior 
has been observed for low fractional fillings typically as- 
sociated with the formation of a static Wigner crystal. 
In practice, however, impurities and defects may influ- 
ence the properties of the rotating crystal (and its exci- 
tations), depending on the magnitude of the excitation 
gap discussed in the previous section. Thus one of the 
main challenges for FQHE observation at such low fillings 
relates to fabrication of high mobility (nearly impurity- 



free) samples 



V. SUMMARY 

In summary, we have carried out the first systematic 
investigations (for 6 < A^ < 9) of structural properties of 
cusp states in parabolic quantum dots at high magnetic 
fields. Our anisotropic conditional probability distribu- 
tions from exact diagonalization show that these states 
are crystalline in character for both low and high frac- 
tional fillings, unlike the liquid-like Jastrow-Laughlin^iiS 
wave functions, but in remarkable agreement with the re- 
cently proposed rotating- Wigner-molecule^ ones. The 
cusp states of A^-electron parabolic QD's are precursors 
to the extended fractional quantum Hall states (and not 
to the static Wigner crystal) due to stabilization of the 
rotating Wigner molecule (having a good angular mo- 
mentum) relative to the static one (that exhibits broken 
symmetry). The rotating Wigner molecule in high B 
does not exhibit global rigidity; instead, it possesses az- 
imuthal rigidity (i.e., all electrons on a given ring rotate 
coherently), with the rings, however, rotating indepen- 
dently of each other. 

Furthermore, we demonstrated pertinent features of 
the spectrum of quantum dots in high B, showing that 
both the ground state and the first excited state corre- 
spond to magic angular momenta (cusp states). For a 
given B, this leads to the appearance of a special gap 
that is not a mere consequence of the finite size of the 
system (and thus it is expected to maintain in the ther- 
modynamic limit, underlying the incompressibility of the 
electron system). Finally, we discussed in detail issues 
pertaining to the implications of the rotating- Wigner- 
molecule theory for FQHE systems at the thermody- 
namic limit (see sect. IV). 

We thank M. Pustilnik for comments on the 
manuscript. This research is supported by the U.S. 
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